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Abstract 

We obtain a general parametrization for the neutrino Major ana mass matrix in which all possible 
deviations from tribimaximal mixing are given by three complex parameters, and the deviation from 
each tribimaximal mixing eigenvector is identified with precisely two of them. This parametrization 
provides a useful tool for classifying the corrections to exact tribimaximal mixing in flavor models. 
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Introduction. Neutrino mass and mixing data to date [1!] strongly favor a "tribimaximal" 
(TB) mixing Q relationship between the neutrino mass and flavor eigenstates. The unitary 
TB neutrino mixing matrix is given by 
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assuming a conventional choice ofphases. Utb is the special case of the general (Pontecorvo)- 
Maki-Nakagawa-Sakata (MNS) [3J] form 
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where s 
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sin 9 



cosOij, for which s± 3 = (and hence 5 is irrelevant), s\ 2 = l/\/3 
and C12 = a/2/3, s 2 3 = c 2 3 = l/y/2, and the Majorana phases /3i j2 , which do not contribute 
to neutrino oscillation processes, are set to zero. While deviations from TB mixing can be 
described in terms of the angles of the MNS matrix, it is more convenient to work directly 
with the eigenvectors of the neutrino mass matrix (i.e., the columns of Utb), which do not 
depend on the conventions for parametrizing unitary matrices. 

The purpose of this note is simply to show that deviations from the tribimaximal form 
of the neutrino Majorana mass matrix are given uniquely by a set of three independent 
complex parameters with the property that setting any two of them to zero preserves one of 
the vectors forming the columns of Utb- Moreover, this property remains true even when 
the nonzero parameter is arbitrarily large. These special parameters are easily identified by 
the specific form of the neutrino mass matrix in any given model. In the physical case, for 
which all the deviations from TB mixing are small, the conventional angles and phases of 
the MNS matrix all have simple forms at linear order in these parameters. We compare our 
approach to other relevant parametrizations at the end of the discussion. 

Diagonalization. The MNS matrix brings the neutrino Majorana mass matrix Mll to 
normal form, diag(mim 2 m 3 ) = U^ NS MllU mns , with m ; > 0. In general, Mil = M is a 
complex symmetric matrix, meaning that its "diagonalization" by means of eigenvectors is 
not completed by the usual means (as evidenced by the presence of C/^ NS rather than ^m NS ), 



but the solution procedure is well known |4j: One solves 
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Mvi = rriiV* , 



(3) 



where m« may always be chosen real and nonnegative, making the normalized vectors Vi 
unique up to overall signs. In the nondegenerate case, one may obtain an equivalent solution 
by multiplying Eq. ([3]) on the left by M*\ 



M*Mvi=m;v,; . 



(4) 



Here the Vi, which are now true normalized eigenvectors of M*M, are defined only up to 
overall phases; nevertheless, since the neutrino mixing angles and Dirac phase 5 are rephasing 
invariant, solving Eq. (J3J) produces all observables except the Majorana phases. Noting that 
M* = M< for complex symmetric matrices, one sees that {mj} are simply the modulus 
squares of the eigenvalues of M. 

While the previous discussion is moot in the special case of exact TB mixing (since U^b is 
real), it is relevant for studying models that deviate from the TB limit. Using the explicit TB 
form Eq. (|Tj) , it is easy to check that all diagonal matrices are in a one-to-one correspondence 
with mass matrices M of the form 
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with complex eigenvalues x — y, x + 2y, and x — y + 2v [6]. This form has three independent 
complex degrees of freedom, x, y, and v, whereas the most general 3x3 complex symmetric 
matrix, has six parameters. The masses m lj2 3 in normal form are the absolute values of 
these eigenvalues; M in this form produces Majorana phases in Umns, as discussed below. 

The unique parametrization for the three remaining parameters, which we call e, Si, and 
82, such that introducing each one of them in turn and keeping the other two zero precisely 
preserves one of the vectors forming the columns of Utb, is given by 
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The six parameters are given by 
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+-(M 12 + 2M 13 + M 22 - M 33 ) , 
--(6M n + M12 - M 13 - 5M 22 - M 33 ) , 



+-(2M 12 - 2M 13 - M 22 + M 33 ) , 
-i(M 12 -M 13 + M 22 -M 33 ), 
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This starting point, a fully general parametrization of Mll in terms of all possible deviations 
from TB mixing, is also the starting point in Ref . , which discusses the preservation of TB 
eigenvectors. Here, however, the preservation of TB eigenvectors is used as an organizing 
principle in defining the parametrization. It is straightforward to check that the unique 
choice of parameters Si, e, and <5 2 preserves the vectors given by the first, second, and third 
columns of £7tb in Eq. (pQ), respectively. In other words, e ^ but Si = S 2 = in this 
parametrization preserves trimaximality, while Si, <5 2 parametrize all possible departures 
from it, and # 2 ^ but e = Si = (/i-r symmetry) preserves bimaximality, while e and Si 
parametrize all possible departures from it. Introducing Si 7^ but e = <5 2 = preserves the 
remaining orthogonal direction (2, —1, — 1) T . 

Parametrizing M as in Eq. ([6]) rather than Utb is more convenient in the construction 
of flavor models, which typically provide textures of elements for Yukawa and Majorana 
matrices rather than the matrices that diagonalize them. Of course, M here refers to the 
Majorana mass matrix Mll in the mass basis for the charged leptons (effectively, taking Yl 
diagonal). In models that specify nondiagonal forms for both Yl and Mll, one faces the 
tricky question of why the mismatch between the matrices that diagonalize them nevertheless 
produces the TB angles to a high degree of accuracy. 

The form of M in Eq. ([6]) also allows for a straightforward calculation of the neutrino 
masses and mixing angles in a perturbative expansion simultaneously linear in e, Si, and 5 2 , 
which is appropriate due to the closeness of the physical values to the TB solution. Solving 
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Eq. (J4]) for the given M produces the results 
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and the corresponding modified eigenvectors 
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Re[y(2x + y)*} ' ^ 

In models in which all the parameters are real, the combinations in Eq. (ITOl) simplify to 

e = e/4t>, 5i = —Si/ (3y — 2v), and 5 = 5 2 /9y. In this case, the independence of the eigenvectors 

on x follows simply from the fact that it enters M as a multiple of the identity matrix. 

These expressions also elucidate the precise conditions for small departures from the TB 

form: |e|, \$i\, \5 2 \ 1. One can easily check that the pathological cases, e.g., v = but 

e^O, lead to degenerate masses, and Um~ns is no longer guaranteed to assume a TB form at 

leading order. 

Now, taking into account the overall phase of appearing in the solution of Eq. ([3]) but 
not Eq. pi), one assembles the new MNS matrix 
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where P = diag{exp(— ~Arg//j)}, with //j being the complex mass parameters in Eq. (jSj): 
mi = ?7mns reduces to TB form with nontrivial Majorana phases when e, Si, and 82 
are set to zero; however, one also notes a different phase convention from Eq. (j2J), in which 
C^mns,23 and £7mns,33 are manifestly real. Multiplication of Umns on the left by an arbitrary 
diagonal phase matrix {i.e., rephasing its rows Umns, %j ~ y ^{/mns, ij), which changes none 
of the observables, resolves this discrepancy. To first order in e, Si, and 82 , one obtains the 
Euler angles 

S13 = \p2\5i - e| , 

S23 = -^|1 -e- 25x| ~-^=[l-Re(e + 250], 

S12 = -^=|1 - 25 2 | ~ -±={1 - 2ReS 2 ) . (12) 

As expected, if bimaximality is preserved (e = 5i = 0), then Sx3 = 0, and the atmospheric angle 
value is given by S23 = 1/ v2 exactly. One also notes that 82 alone in this approximation 
parametrizes departures from the solar angle value Sx2 = V y/3. Using global fit values of the 
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one finds that the parameters are indeed small: 



\Si - e\ = 0.07+ 006 



0.04 5 



Re(e + 25!) = 0.00±g;g? , 

Re5 2 = 0.023±g;gg. (14) 

The rephasing that leads to real Umns, 23 and Umns, 33 is 

Ot = ^Arg(^ 3 ) + m , where 7^,2,3 = 0, +Im(e + 2^), -Im(e + 2^) . (15) 

The MNS phases may then be obtained from Eqs. (!9l)- (jTTl) : 

5 = -Arg(5x - e) , 

Pi = g ( Ar S^ 3 ~ Ar S^i) _ Im ^2 , 

P2 = ^(Arg / u 3 -Arg/i 2 )-2Im5 2 , (16) 



and the Jarlskog parameter is given by 

J=ilm(5 1 -e). (17) 

Note that (3i, @2, and 8 can all be nonzero even when e, 5i, 5 2 — > 0. Although the mass-squared 
splitting ratio Am^/Am^ = 0.0319^oooi8 is experimentally small Q], it is nonvanishing in 
the TB limit: 

Ami 3Re [y(2x + y)*} ~ ~ 

Am| 2 " Re [(3y - 2v)(2x + y + 2v)*] + 1 ' u 2) ' { } 

Discussion. A number of parametrizations of the MNS matrix have appeared in the 
literature [8|. Two recent parametrizations in particular describe deviations from TB in 
terms of three (real) parameters plus the CP-violating phase 5. The relations of these two 
sets, r, s, and a in Ref. [sj and ei 3 , e 2 i, and e 32 in Ref. [h]], to those defined here at linear 
order are given by 

\x ~\ r 613 
Re (e + 2Si) = -a = -e 32 , 

^4 = -£ = -*. (is) 

The advantages of the current parametrization have already been identified, but to sum- 
marize: (1) The parameters e, Si, 82 have a natural interpretation in terms of which TB 
eigenvector each preserves, even when the chosen parameter is not perturbatively small. (2) 
The mixing angles and Dirac phase of the MNS matrix have simple expressions at linear 
order in e, 8\, <5 2 . (3) Our parameters allow a direct analysis of mass matrix corrections that 
lead to a deviation from exact TB mixing in flavor models. For example, in one well-known 
model [Hi], the entries of Mn are determined by four vacuum expectation values (yevs) of 



Higgs triplets, labeled a, b, c, d. From the mass matrix given in Eq. (19) of Ref. [11] . it is 
straightforward to check [using, for example, Eq. ([7D] that 

8 X = 6 2 = and e = b - c . (20) 

Our parametrization shows immediately that this model is always exactly trimaximal. The 
deviation from bimaximality is reflected in s 13 = 7^ | | when \b — c| |c — a\ and all vevs 
are real. 
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